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■ Abstract. FUSE far-UV spectroscopy of the nearby metal-deficient {Zq/8) 

cometary Blue Compact Dwarf (BCD) galaxy Mrk59 have been used to investigate 
element abundances in its interstellar medium. The lack of diffuse H2 (A'^(H2)< 
lO^^cm"'^) is due to the combined effect of a strong UV radiation field and a low 
' metallicity. By fitting the line profiles with multiple components having b = 7 km 

^SJ ' s~^ and spanning a radial velocity range of 100 km s~^, we derive heavy element 

abundances in the neutral gas. We find log N(0 l)/N{Y{ l) = - 5.0±0.3, about a 
factor of 10 below the oxygen abundance of the supergiant H 11 region. 
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cn . 1 FUSE Observations of Mkr 59 

m 

0\ ; The Blue Compact Dwarf (BCD) galaxy Markarian 59 (Mrk 59) ee I Zw 49 

' belongs to the class of cometary BCDs with an intense starburst at the end 

of an elongated low surface brightness stellar body. This BCD have numerous 
knots composed of a chain of H ii regions, resulting from propagating star 
formation along the galaxy's elongated body and ending with the high-surface- 
brightness supergiant H ii region. From deep ground-based spectrophotometric 
O \ observations of the supergiant H ii regions Noeske et al. (2000) derived an 

' Oxygen abundance log 0/H — -4.011 ± 0.003 {Zq/%). O abundances were 

, also derived for two other emission knots along the elongated body and were 

found to be the same within the errors. The small scatter in metallicity along 
the major axis of Mrk 59 (~ 0.2 dex) suggests that the mixing of elements in 
' the ionized gas has been efficient on a spatial scale of several kpc. 

^ \ Mrk 59 was observed during 7,865 s on 2000, January 11 with FUSE (Moos 

et al. 2000). The LWRS large entrance aperture (30"x30") was used, so that 
all of Mrk 59 which is ~ 20" across (Noeske et al. 2000) is included within it. 

Mrk 59 shows numerous interstellar absorption lines from the H i Lyman 
series (Ly /3 to Ly A (H i 11)) and from other atoms and ions such as C ii, 
C III, N I, N II, N III, O I, Si II, S III, S IV, Fe ii and Fe iii. Several stellar 
features such as the Si iv and S vi lines are also detected (Thuan et al. 2001a). 



2 Narrow cores and broad wings of the H I Lyman series 
line profiles: interstellar and stellar absorption 

We found that the shapes of the H i Lyman series lines cannot be reproduced 
with Voigt profiles characterized by a single H i column density. Damping 



wings are clearly visible, signaling a very large H i column density ~ 5 x 
10^^ cm~^. However, this high column density is not consistent with the 
shape of the strongest H i lines (Lyman (3 to Lyman S). For instance, although 
damping wings are also present, the core of the Lyman 6 line is thin and cannot 
be fitted with a model profile with the column density needed to fit the higher 
lines in the series. The thin core implies a smaller column density, of the order 
of 10^^ cm~^. We interpret that by a non-interstellar but stellar origin of the 
wide damping wings. We conclude that the broad wings of the above lines 
arise in the photospheres of numerous B stars, while their narrower cores are 
caused by interstellar absorption. By fitting the cores of the Lyman series 
lines, we obtain an interstellar H l column density of ~ 7x10^° cm~^. 

3 Upper limits on the diffuse H2 content of Mrk 59 

No line of H2 is seen at the radial velocity of Mrk 59. Using nine H2 Lyman 

bands (0-0 to 8 0) , we conclude that the total column density of diffuse H2 is 
<10^^cm~^. This implies that the ratio of H2 to Hi is <10~^. The correspond- 
ing average molecular fraction f=2N{a2)/{N{R i)-|-2Ar(H2)) is < 3 x IQ-^^. 

By calculating the amount of H2 molecules formed on the surface of dust 
grains, we find that the expected iV(H2) must be ~2xl0^^ cm~^, consistent 
with our upper limit. As compared to the Milky Way, the low column density 
of diffuse H2 in Mrk 59 is due to the combined effects of a large UV flux which 
destroys the H2 molecules and of a low metallicity resulting in a scarcity of 
dust grains on which H2 form (see the full discussion in Thuan et al. 2001a). 

4 Heavy element abundances 

To estimate the abundances, we consider two cases. The first case presented 
during the meeting assumes that there is a single velocity component along 
the thousands of lines of sight to the stars of Mkr 59. In this case, absorption 
lines that are broader than the instrumental Point Spread Function and do 
not go down to zero arc supposed to be not saturated. The second case, 
suggested during the meeting (York, public communication), assume multiple 
velocity components along the multiple lines of sight, some of which may have 
saturated profiles. 

4.1 Profile fitting with a single velocity component 

We consider first the single interstellar velocity component case. Some of the 
lines, like the Fein and Sin lines, are very broad with unreasonably large 
Doppler widths (6>100 km s"-'^), suggesting, as for the Hi lines, that these 
metal lines have some contamination by stellar absorption. The resulting col- 
umn densities (log A^i) and Doppler widths (61) along with their error bars are 
given in Table 1. The error bars of the column densities include the uncertainty 



Table 1: Heavy element column densities in Mrk 59 
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N I 
I 
S IV 
Fe III 


40±20 

40±20 
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a) Fitting with a single Voigt profile, bi is in kms~^. Error bars are 3 a limits. 

b) Fitting with multiple velocity components, each with 62 = 7 km s~^, and 
spanning a velocity range of 100 kms~^. Error bars are 2 a limits. 



in b. These abundances of the metals relative to A'^(Hi) arc extremely low. In 
particular, with log iV(0 i)/iV(Hi) = - 5.7, the oxygen abundance in the Hi 
gas could be ~ 50 times lower than the oxygen abundance in the ionized gas, 
as determined from the optical emission line spectrmn of the supergiant Hll 
region in Mrk 59. However, these very low abundances may be underestimates 
due to the saturation of some lines of sight. 

4.2 Profile fitting with multiple velocity components 

We are observing through the FUSE aperture along thousands of lines of sight 
towards UV-bright stars. The observed spectrum is the sum of many spectra. 
It can happen that some lines of sight have saturated absorption lines with 
a small b parameter, but because they have different radial velocities spread 
over several tens of kms^^, the broad absorption line resulting from the sum 
of many narrow absorption lines does not go to zero intensity, and its width is 
larger than the Point Spread Function. In that single velocity compo- 

nent fit to the line profile would result in an overestimate of the 6 parameter 
and in an underestimate of the column density. 

To investigate this issue, we have calculated profiles resiilting from thc^ 
addition of multiple lines of sight. We adopt the simple model where the 
different lines of sight have radial velocities distributed uniformly between 
^'Mrksg — Av/2 and iiMrk59 + Av/2, where Av is the spread in velocity due 
mostly to the velocity dispersion of the system of absorbing clouds along the 
multiple lines of sight, and also partly to the wavelength smearing caused by 
the extension of Mrk 59 within the aperture. We determine the best estimates 
of b (assumed to be the same for all lines of sight) and Av by the simultaneous 
fit of several lines of O i, N I and Fe 11. Each of these elements having several 
transitions with different oscillator strengths, this well constrains b and At; 
to reasonable values: 6=7^3^ km s~^ and At;=100±20kms~^ (2cr limits, see 
Fig 1). Note that the value of Av is very close to the FWHM of 92kms~-^ 
of the H I profile (Thuan et al. 2001b). The heavy element column densities 
derived with the multi-component fit are given in Tablet (log7V2). Those are 
larger than in the case of a single component fit, by a mean factor of '--^5. 
The very large column densities derived for C III, S III and Fe iii are caused by 




Figure 1: Fit to the O i 1039 A absorption line with a column density iV(0 l) 
= 7 X lO^^cm"^. The fit is obtained by the addition of multiple profiles. Av 
and b parameters are constrained by the simultaneous fit of several lines of O i, 
N I and Fe ii. The thick solid line gives the resulting fit to the data with Av = 
100 km s~^ and 6 = 7 km s^^. The dotted line represents an example of one 
line of sight with 6 = 7 km s~^. With this column density, such an individual 
line of sight is saturated and goes down to the zero intensity level. However 
the absorption line resulting from the addition of many narrower profiles does 
not go to the zero intensity level, although it is broader than the width of the 
instrumental Point Spread Function (~ O.lA). 

stellar contamination. It is interesting to note that, even with a higher derived 
oxygen abundance, log A'^(0 i)/7V(H i) = - 5.0±0.3, the H i absorbing cloud 
still has a metallicity lower than that of the supergiant H ii region by a factor 
of ~ 10. This suggests self-contamination of the H ii region by heavy elements 
released during the present burst of star formation (Kunth & Sargent 1986). 
While mixing of these newly formed heavy elements appears to have occurred 
on the scale 2 kpc) of the H ii regions (Noeske et al. 2000), it has not 
had time to occur for the whole neutral gas component as the H i envelope 
surrounding the star-forming regions is much more extended (its diameter is 
~ 19 kpc from the VLA map by Thuan et al. 2001b). 
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